The birth place of the type Ic Supernova 2007gr 
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ABSTRACT 

We report our attempts to locate the progenitor of the pecuhar type Ic 
SN 2007gr in HST pre-explosion images of the host galaxy, NGC 1058. Aligning 
adaptive optics Altair/NIRI imaging of SN 2007gr from the Gemini (North) Tele- 
scope with the pre-explosion HST WFPC2 images, we identify the SN position 
on the HST frames with an accuracy of 20mas. Although nothing is detected at 
the SN position we show that it lies on the edge of a bright source, 134 ± 23mas 
(6.9 pc) from its nominal centre. Based on its luminosity we suggest that this 
object is possibly an unresolved, compact and coeval cluster and that the SN 
progenitor was a cluster member, although we note that model profile fitting 
favours a single bright star. We find two solutions for the age of this assumed 
cluster; 7 =F O.SMyrs and 20 — SOMyrs, with turn-off masses of 28 ± 4M0 and 
12 — 9Mq respectively. Pre-explosion ground-based i^'-band images marginally 
favour the younger cluster age/higher turn-off mass. Assuming the SN progenitor 
was a cluster member, the turn-off mass provides the best estimate for its initial 
mass. More detailed observations, after the SN has faded, should determine if 
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the progenitor was indeed part of a cluster, and if so allow an age estimate to 
within ~2 Myrs thereby favouring either a high mass single star or lower mass 
interacting binary progenitor. 

Subject headings: stars: evolution - supernovae:general - supernovae: individ- 
ual(2007gr) - galaxies: individual(NGC 1058) 



1. INTRODUCTION 



Direct detection of the progenitors of core-collapse supernovae (CCSNe) in archival 
observations of nearby galaxies, particularly those from the Hubble Space Telescope (HST), 
has recently proved successful. Luminous red supergiant progenitors of the hydr o gen rich 



type II-P SNe have now been identified (jSmartt et al.l l2004l : iMaund et al.l 12005 



Li et al. 



20061 ). with initial masses above 8Mq. However the progenitors of type Ib/c SN, believed to 
be t he cores of massive stars that have lost their outer envelopes, have so far eluded discovery 
(e.g. iMaund fc Smarttll2005l ). The two nearest type Ib/c SNe with deep pre-explosion images 
('SN2002a p in M74 at ~9Mp c and SN2004gt in NGC4038 at ~14 Mpc) had no progenitor 



detected (IMaund et all l2005l : iGal-Yam et all l2005l : ICrockett et all 120071 ). The progenitors 
were limited to either lower mass stars that were stripped of their hydrogen-rich envelope 
via binary interaction, or single stars with masses > 30 — 4OM0 for which radiatively driven 
mass loss may have been sufficient to remove the outer layers. In the case of the Type 
lb SN 2006jc the progenitor was discovered in a n LBV-like outburst s tate 2 years before 
explosion, suggesting it was a very massive star (jPastorello et al.l 120071 ). Some fraction of 
Type Ib/c SNe also produce g amma-ray bursts, hence any information on their progenitors 
would be critical for the field (IWoosley fc Bloom! |2006| ). 



In this letter w e study the explosion site of SN 2007gr in NGC 1058 (1 0.6 ± l.SMpc: 
Pilvugin et al.ll2004l ) and confirm it as a type Ic. SN 2007gr was discovered bv lMadison fc Li 



(2007) on 2007 Aug 15.15 UT, located at ajaooo = 2M3"^27^98, 5j2ooo = +37°20'44".7. 
( 120071 ) spectroscopically classified SN 2007gr as Type Ib/c, but suggested that the presence 
of strong O I A7774 and Ca II IR triplet absorptions were similar to those features observed in 
Type Ic SNe. We adopt a total reddening towards SN 2007gr of E{B -V) = 0.08 (Valenti et 
al. in prep). Where standard filter sets were employed extinctions were calculated using the 
laws of Car delli, Clayton fc Mathi s (1989), and for HST observations using the Ax/ E{B — V) 
relations of I Van Dyk et al.l (Il999l ). 



Chornock et al, 
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2. OBSERVATIONS AND DATA ANALYSIS 

The pre-explosion site of SN 2007gr was observed on 2001 Jul 03 with the Wide Field and 
Planetary Camera 2 (WFPC2) on board HST as part of programme GO-9042. Observations 
were made in the F450W and F814W filters (460s e ach) and phot ometry was carried out 



using the PSF-fitting photometry package HSTphot flDolphinll2000l ). Further ground-based 
observations of NGC1058 were found in the 4.2-m William Herschel Telescope (WHT) and 
the 2.5m Isaac Newton Telescope (INT) archives. The WHT INGRID fC— band images from 
2001 Oct 6 had a total exposure time of 1320s and an image resolution of 1'.'2. The INT 
WFC images from 2005 Jan 13 were taken using r'-band (180s) and narrow Hq, filters (600s) 
with image quality of 0'.'9. 

Ground-based adaptive optics (AO) observations of SN 2007gr were taken on 2007 Aug 
19 using Altair/NIRI on the 8.1-m Gemini (North) Telescope. Exposures totalling 930s 
were acquired in i^'-band using the SN as the natural guide star. The reduced image is of 
exceptionally high quality showing near diffraction limited resolution of 0'.'08, with a pixel 
scale of 0'.'022 providing excellent sampling of the PSF. Fig. [1] shows that this Gemini image 
even resolves some objects which appear blended in the HST data, an effect of the O'.'l pixel 
scale of the WFPC2 chip. In order to confirm the classification of the SN, a spectrum was 
obtained on 2007 Sept 11 using the WHT equipped with ISIS (R300B+R158R). 

The site of SN 2007gr fell on the WF2 chip of the pre-explosion WFPC2 observations. 
To determine the position of the SN site accurately on these images we calculated a trans- 
formation to map the coordinate system of the post-explosion Gemini image to the WFPC2 
pre-explosion frames using the IRAF geomap task. Positions of 19 stars common to both 
the WFPC2 F814W and Gemini i^'-band images provided a general geometric transforma- 
tion (involving shifts, scales and rotation) with an RMS error of ±20mas. The SN position 
was measured in the Gemini image and transformed to the coordinate system of the F814W 
pre-explosion frame yielding a pixel position on the WF2 chip of [250.65, 118.34]. Fig.[T^ & 
b show that the SN site is very close to, but not coincident with, an object that is bright 
in both the F450W and F814W images. Its position was measured as [250.10, 119.56] with 
an accu racy of ±llma s. The SN site is therefore 0'.'121 west and 0'.'056 south of this bright 



object. iLi et al.l (120071 ) give a position for SN2007gr on the WFPC2 images that is around 
66 mas further north than our position, meaning our respective positions are not coinci- 
dent within the combined errors (combined error from this study ±23mas, Li et al. claim 
±20mas). We suggest that the Li et al. error has been underestimated given that they 
used a MEGACAM image of the SN from the Canada France Hawaii Telescope to perform 
differential astrometry, which is of much lower resolution than our Gemini data. 



The position of the SN on the INT and WHT pre-explosion observations was determined 
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using a similar method. The INT r'-band image (Fig.[T]i) shows significant flux at the SN 
position, but this is almost certainly due entirely to the resolved sources observed in the 
HST frames. Subtraction of the r'- band continuum from the image identifies regions of 
strong Ha emission. Since hot OB-type stars are required to ionise the hydrogen gas, 
emission is generally associated with clusters of young, massive stars. The image subtraction 
package ISIS 2.2 (Alard & Lupton 1998; Alard 2000) was used to match the PSF, fiux and 
background of the INT r' image to its counterpart. Fig.[T^ shows several regions of quite 
strong emission surrounding the SN site, but none exactly coincident with its position. In 
the WHT K-band image (Fig.[T]F) nothing is visible at or near the SN position. 



OBSERVATIONAL RESULTS AND DISCUSSION 



SN2007gr occurred in a large, extended star-forming region that in Fig{T] is at least 
5" in diameter (250pc). Within this there are bright sources, some of which are point- 
like and others which are clearly extended. Most of these are likely to be compact stellar 
clusters rather than single individual supergiants. Compact star clusters tend to be clustered 
together into large star forming complexes of the order 50-300 pc, and the recent study of 
the cluster comple xes in M51 shows m any that are morphologically similar to the region 
hosting SN2007gr (IBastian et al.ll2005l ). HSTphot photometry of the pre-explosion WFPC2 
observations gave apparent magnitudes for the bright object close to the SN site of F450W 
= 21.47±0.02 and F814W = 20.89±0.02. Converting these to absolute standard magnitudes 
using the dist ance to NGC 105 8 and extinction corrections (Section [1]), and the colour 
corrections of iMaund fc SmarttI J2005h . gave Mb = -8.99, and (B - I) = 0.42. If it is 
a single star it is a 3O-4OM0 post main sequence supergiant star, or a more massive but 
extinguished {Ay ^ 1 mag) main seque nce object. Such bri ght, single stars are not without 



precedent in star-forming galaxies (e.g. iBresolin et al.ll200ll ). but they are few in number in 



any typical spiral due to the short lifetimes of the blue- yellow supergiant phase (e.g. Eldridge 
& Tout 2004). 

Surrounding the SN location are several other objects of similarly high luminosity that 
have colours consistent with O- to F-type supergiants (Fig. [2]). It is unlikely that one would 
observe such a large number of extremely massive stars d uring this very short evolutionary 
supergiant phase in such a small surface area of the disk. iBastian et al.l (120051 ) suggest that 
sources brighter than My < —8.6 are most probably clusters. Several of the objects appear 
extended implying that they are most probably not single stars. Most however, including 
the bright source nearest the SN location, are not significantly broader than the stellar PSF 
(FWHM~1.5 pixels (0'.'15)), which at a distance of 10.6 Mpc corresponds to 7.7pc. Given 
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its absolute magnitudes and its proximity to the SN site (~6.9pc) one could suggest that 
this bright object is a coeval cluster of diameter less than ~7.7pc, which hosted SN 2007gr. 
If it could be confirmed as such a cluster this is potentially a powerful way of determining 
the initial mass of the SN progenitor. 

In an attempt to constrain the true spatial extent of this object and others we used the 
Ishape routine, part of the BAOlab package (Larsen 1999, 2004). Larsen (2004) suggests 
that Ishape can limit the intrinsic size of a significantly detected source (S/N ~50cr) down to 
~10 per cent of the PSF FWHM, which here corresponds to ~1 pc. Elliptical Moffat profiles 
with a power law index of -1.5 were convolved with TINYTIM (Krist & Hook 1997) WFPC2 
PSFs and fitted to all the objects within 3" of the SN position. Those sources which are 
detected at >10cr level and have fitted FWHM values >0.15 pixels (10 per cent of the PSF 
FWHM) in one or both of the HST observations are identified as extended in Fig. [21 The 
bright source nearest the SN position is best fit with a delta function, suggesting that its 
diameter is no larger than ~1 pc, and that it is more likely a star rather than a cluster. 

However, aperture photometry of the Gemini AO image, calibrated using the two bright 
2MASS sources north and south of the SN, revealed that several sources around the SN 
location have significant K-band excesses when compared to main sequence and supergiant 
stellar colours (see Fig.[2]), implying that they are not individual massive stars. Two of these 
sources are measured as point-like by Ishape. Unfortunately we cannot measure a K-band 
magnitude for the object nearest the progenitor site in the Gemini image due to the brilliance 
of the SN. Additionally the Ha image (Fig.[T^) gives some idea of the location of the youngest 
OB-stars within this large star-forming region. Sources in the HST images coincident with 
the Ha flux are marked in Fig.[2l Their association with significant Ha emission implies that 
they are most likely young compact star clusters rather than single massive supergiants. We 
note, however, that no significant Ha emission is observed exactly coincident with the object 
nearest the SN site. 

If this object is actually a massive single star rather than a cluster, and the progenitor is 
an unrelated star within this extended star-forming region, then we can attempt to derive a 
luminosity limit for the undetected SN progenitor. We added fake stars of various magnitudes 
into the WFPC2 images at the position of the SN progenitor and tried to recover both the 
real and fake sources simultaneously with PSF-fitting techniques. We defined the detection 
limit in each filter as the magnitude of the fake source that must be added to the image in 
order for it to be independently detected at a position coincident with the progenitor site. 
The detection limits were thereby found to be F450W > 23.7 (Mf45ow > -6.7) and F814W 
> 21.7 (Mf814w > —8.6). The detection limit for the WHT K-band image was found to be 
K > 19.7 {Mk > —10.45) by adding fake stars of increasing luminosity at the SN position 
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and determining, by eye, at what magnitude an obvious source became visible. 

Fig.O shows our WHT + ISIS spectrum compared with those of well studied type Ib/c 
SNe. The phase corresponds to about two weeks past maximum (Valenti et al. in prep.). 
While in SN 1999dn the He I lines are prominent, they are weak (if at all visible) in the type 
Ic SNe and also in SN 2007gr, suggesting that SN 2007gr is of type Ic. SN 2007gr shows 
clear similarity to SN 19941 (and, more marginally, to SN 2004aw). However, the presence 
of narrower features due to metals in the spectrum of SN 2007gr makes this object rather 
peculiar. 



4. IMPLICATIONS FOR THE PROGENITOR OF SN2007gr 



If we assume that the bright object closest to the SN site is a single, massive star, 
one can attempt to place constraints on the properties of the SN progenitor using the de- 
tection limits derived in the previous section. However, these r nagnitude limits are much 
shallower t han those found fo r the progenitors of SN 2002ap (jCrockett et al.l 120071 ) and 
SN 2004gt flMaund et al.ll2005l ). Since the SN was of Type Ic we can reasonably assume that 
its progenitor was a Wolf-Rayet (W-R) star (e.g. Eldridge & Tout 2004). Comparing our 



m agnitude limit o : 



Me450W > —6.7 with B-band photometry of W-R stars shown in Fig. 8 



of ICrockett et al.l (120071 ). we see that almost all of these observed W-R stars fall below our 
detection limit. 

If, however, the bright object closest to the SN site is an unresolved stellar cluster, then 
it must have a diameter of less than ~7.7pc (FWHM of the PSF) or perhaps less than ~1 pc 
{Ishape analysis). In the case of the latter the progenitor would then be on the outskirts of the 
cluster, although still quite possibly coeval, e.g. NGC 3603 in the Milky Way. If SN 2007gr 
arose from one of the most evolved stars in such a coeval cluster then the cluster age would 
give an estimate of the age and initial mass of the progenitor. A similar argument has been 



used to estimate t he progenitor mass of the type II-P SN 2004dj (iMaiz-Apellaniz et al.ll2004 



Wang et al.l l2005l ). through starburst99 p opulation fitting of t he UBVRIJHK spectral 
energy distribution. The starburst99 code (ILeitherer et al.lll999l ) computes a coeval stellar 
population at a range of ages, with a user-defined initial mass function. A cluster of mass 
~ 4000 Mq at an age of 7 ^ 0.5 Myrs produces the absolute Mb = -8.99 and 5 - / = 0.42 
(assuming solar metallicity, an upper mass limit of IOOMq, a Salpeter IMF, and a coeval 
burst of star-formation). The Geneva stellar models imply a turn-off mass of 25M0 for the 
cluster, while our STARS stellar models (Eldridge & Tout 2004) imply a similar mass of 
28 ± 4M0. These are compatible with th e minimum stellar mass thought to produc e W-R 



stars through single stellar evolution (e.g. iMassey et al.l 1200 ll : lEldridge &: Toutll2004l ) 
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However this age solution is not unique as starburst99 can produce B — I^ 0.4 when the 
cluster has evolved to 20—30 Myrs. At this age the turn-off mass would be between 12 — 9Mq, 
much too low to produce a W-R star via single stellar evolution, and we would be forced 
to conclude that the progenitor was a lower mass W-R star resulting from an interacting 
binary. This degeneracy can only be broken with the extra wavelength coverage of the host 
cluster. In particular the ?7-band is an excellent discrirninator of age, and would allow its 
determination to within a few Myrs (IBastian et al.ll2005l : iMaiz-Apellaniz et al.ll2004j ). At the 
older age the cluster becomes much redder in the NIR bands, with an estimated I — K = 1.1 
and hence we would expect to detect the source at an apparent magnitude of K = 19.7. The 
detection limit of the WHT K-band image lies just at this level, and this would favour the 
younger solution, suggesting the W-R progenitor formed from an initially 25-3OM0 star. 

In approximately two years time, the SN will have faded enough that the putative clus- 
ter can be studied in detail, in the UV, optical and NIR with ground-based AO imaging and 
hopefully a fully refurbished HST. With these higher spatial resolution and multi- wavelength 
images of the immediate environment we will be able to determine if indeed it is a com- 
pact cluster and better ascertain its properties. Given the difficulties i n pinning down the 



progenitors of Ib/c SNe, even with very deep pre-explosion imaging (IMaund et al.l 12005 



Crockett et al.l 120071 ). this particular case may offer us the best opportunity to estimate the 
initial mass of a type Ic progenitor. The proximity and luminosity of SN 2007gr will make 
it one of the best studied, and modelled type Ic SNe to date, and it will be interesting to 
compare future ejecta mass estimates with our cluster turn-off mass. 



Observations were from : the NASA/ESA Hubble Space Telescope obtained from the 
Data Archive at the Space Telescope Science Institute; the Gemini Observatory; the Isaac 
Newton Group Archive from CASU Astronomical Data Centre. We acknowledge funding 
through EURYI scheme, ESF (SJS, RMC), NSF/NASA grants AST-0406740/NNG04GL00 
(JRM), Academy of Finland project: 8120503 (SM). 
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Fig. 1.— Pre- and Post-explosion images of the site of SN 2007gr in NGC 1058. Each of the 
images is centred on the SN position and orientated such that North is up and East is to 
the left. The position of SN 2007gr is located at the centre of each of the panels indicated 
by the cross hairs, a) Pre-explosion HST + WFPC2 F450W image, b) Pre-explosion HST 
+ WFPC2 F814W image, c) Post-explosion Gemini + Altair/NIRI K-band image with the 
SN clearly visible, d) Pre-explosion INT H- WFC r'-band image, e) Pre-explosion INT -|- 
WFC Ha minus r' image, f ) Pre-explosion WHT -|- INGRID K-band image. 
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Fig. 2. — Colour-magnitude diagram of sources within a radius of approximately 3" (150 pc) 
of the SN position. Photometry is from the WFPC2 F450W and F814W pre-explosion 
observations. All sources are marked by a filled circle. Objects which Ishape measured as 
extended are marked with triangles, those with K band excesses of >0.5 mags when compared 
to single stellar colours are marked with crosses, and those associated with strong Ha emission 
are indicate d by open circl e s. Th e bright object closest to the SN position is marked with 
a filled star. iBastian et al.l (120051 ) suggest that anything brighter than My < —8.6 in these 
starforming complexes is likely to be an unresolved cluster rather than a single, massive star. 
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Fig. 3. — Comparison of spectra of the classical Type lb SN 1999dn (Branch et al. 2002), the 
peculiar Type Ic SN 2004aw (Taubenberger et al. 2006), the classical SN Ic 19941 (Filippenko 
et al. 1995) and SN 2007gr at 12-15 days past maximum. The spectrum of SN 2007gr shows 
several unusually narrow P-Cygni lines. No reddening correction was applied to the spectra, 
while they are reported to the host galaxy rest wavelength. The vertical dashed lines mark 
the positions of the minima of the He I lines in the spectrum of SN 1999dn, blueshifted 
by about 9500 km/s. No prominent He I lines are visible in the spectrum of SN 2007gr, 
suggesting a Type Ic classification. 



